Full counting statistics in a propagating quantum front and random matrix spectra 
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One-dimensional free fermions are studied with emphasis on propagating fronts emerging from 
a step initial condition. The probability distribution of the number of particles at the edge of the 
front is determined exactly. It is found that the full counting statistics coincides with the eigenvalue 
statistics of the edge spectrum of matrices from the Gaussian unitary ensemble. The correspondence 
established between the random matrix eigenvalues and the particle positions yields the order statis- 
tics of the right-most particles in the front and, furthermore, it implies their subdiffusive spreading. 



The theory of quantum noise progressed rapidly dur- 
ing the last decades and it became an important research 
area in the studies of transport in mesoscopic systems 
[H, 0]. One of the key concepts of this field is the full 
counting statistics (FCS), giving the probability distri- 
bution of the charge transmitted through a conductor 
Non-Gaussian fluctuations of the FCS carry impor- 
tant information and the measurement of higher order 
cumulants has recently become accessible in a number 
of experiments on tunnel junctions Q, [Bj], quantum dots 

0] and quantum point contacts The FCS has also 
been proposed as a tool to extract the entanglement en- 
tropy [§] in terms of the measured cumulants Toi- 12\. 

Theoretically, the FCS is best understood for the trans- 
port of non-interacting electrons between conductance 
channels separated by a scatterer At zero tem- 

perature, an even simpler setup is to consider a one- 
dimensional system with perfect transmission, where the 
current is induced by preparing an initial state with two 
parts of the system having unequal particle densities. In 
this case, the complete time evolution of the FCS can 
be followed [13, LL4| and interesting results (e.g. particle 
number fluctuations logarithmically growing with time 
HI . 15| ) can be obtained. 

The above setup with a step-like initial density is also 
remarkable because, due to the density bias, a front 
builds up and penetrates ballistically into the low density 
region 16]. This front displays an intriguing staircase- 
like fine-structure in the density [ItJ which has been also 
observed in the fronts of various quantum spin chains 
and in equivalent fermionic systems H|"21|. For free 
fermions, the fine structure appears to be related to 
the curvature of the single-particle dispersion around the 
Fermi points and thus cannot be accounted for in the 
usual semiclassical picture. Nevertheless, the staircase 
appearance of the front suggests a particle interpretation 
nj and thus the FCS appears as a natural candidate for 
exploring the details of quantum fronts. 

Our aim here is to probe the quantum noise in the 
front region by developing the FCS in a frame comoving 



with the edge of the front. For the free fermion case, our 
main result is the discovery of a connection between the 
FCS at the front edge and the distribution functions of 
the largest eigenvalues of the Gaussian unitary random 
matrix ensemble. This result allows us to characterize 
the front in terms of particles whose statistics, and in 
particular, their order statistics can be derived exactly. 

The system under consideration is an infinite chain of 
free spinlcss fermions described by the Hamiltonian 
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where cj^ is a fermionic creation operator at site m. Ini- 
tially, all the sites with to < are filled while those with 
m > are empty and thus the correlations at time t = 
are given by 



b mn TO, n < 

else 



(2) 



Since H is diagonalized by a Fourier transform with a 
single-particle spectrum uj q = — cos q, the time evolution 
of the Fermi operators c m (t) can be obtained in terms of 
Bessel functions J m if) as 
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The step function in the density ri m (t) spreads out bal- 
listically and for m > is given by 16 1 



n m (t) = (<4(i)c„ 
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with snapshots of the density profile shown on Fig. [T] 
The curves for different times collapse in the bulk by 
choosing the scaling variable m/t, and the emerging pro- 
file can be understood by semiclassical arguments ljj . 
However, there is a nontrivial staircase structure emerg- 
ing around the edge of the front, shown by the shaded 
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areas in Fig. [TJ The size of the edge region was found 
to scale with i 1 / 3 and it was argued that each step con- 
tains one particle 17[ . Since the origin of this behaviour 
cannot be captured by a semiclassical treatment, a more 
detailed understanding will be obtained by looking at 
the FCS which includes all the higher order correlations 
of the particle number. 
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FIG. 1: Density profiles n m (t) at increasing times. The edge 
of the front with the emerging staircase structure is shown 
by the shaded region. The profile for m < is given by 
n m {t) = 1 — ni- m (t) and is not shown. 

The FCS is defined through the generating function 



X (X,t) = (exp(iXN A (t))} 



(5) 



where N^t) — J2m£A c rn(t) c m(t) is the particle number 
operator in subsystem A at time t. The subsystem will 
be chosen to include only the edge of the front and thus 
A itself is time dependent. The generating function can 



be written as a determinant 13j, [14 1 



X (A,t)=det[l + (e !A -l)C(i)] 



(6) 



where C(f) is the reduced correlation matrix at time t 
with matrix elements C mn (t) — (cj„(t)c ra (t)) restricted to 
the subsystem, to, n € A, while 1 is the identity matrix 
on the same interval. For the s tep initial condition one 
obtains the matrix elements as 1161 



t 



2(m — n) 



[J m -l(t)J n (t) - J m (t)Jn-l(t)} (7) 



For convenience, we define C mn (t) by dropping the phase 
factor i"~ m in Eq. ([7]) which corresponds to a simple uni- 
tary transformation of the matrix C(i) and thus leaves 
x(A, t) invariant. Note, that C(t) has exactly the same 
form as the ground-state correlation matrix of a free 
fermion chain with a gradient chemical potential 22 1. 



Thus our results can be directly applied to the static in- 
terface problem as well. 



In order to explore the edge region of the front which 
scales as t 1 ^ 3 , we introduce scaling variables x and y 
through to = t + 2- 1 /H 1 / 3 x and n = t + 2" 1 / 3 i 1 / 3 y. 
We also use the identity J m _i(t) = J m {t) + yJ m (t) to 
rewrite C mn (t) in terms of Bessel functions and their time 
derivatives as 
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2(m — n) 



Jm{t)J n {t) - J m (t)J n (t) +-J m {t)J n {t) (8) 



1 



In the asymptotic regime t — > oo, one has the following 
forms 23] for the Bessel function and its time derivative 



J m (t) « 2 1 ' 3 t- 1 / 3 Ai( I ) 
j m (t) « -2^t-^M\x) 



(9) 



and the other pair of equations is obtained by interchang- 
ing to — in and x — > y. Substituting into Eq. (|5J| one 
arrives at 

C mn {t) « 2 1 / 3 t- 1 / 3 ^(a; ; y)+2- 1 / 3 t-2/3 Ai(a . )Ai(y) (10) 



where 



K(x,y) 



Ai(x)Ai'(y) - M'(x)M(y) 



(11) 



The factor 2 1 / 3 t~ 1 / 3 in the first term of Eq. ([TP} accounts 
for the change of variables while the other term vanishes 
in the scaling limit. Therefore, the reduced correlation 
matrix is turned into an integral operator and the gener- 
ating function is given by the Fredholm determinant 



X (X,s)=det [1 + (e lA - l)K] 



(12) 



where s denotes the scaled left endpoint of the domain 
A = (s, oo) over which the kernel K (x, y) is defined. 

The Airy kernel (fTTj) is well known in the theory of 
random matrices from the Gaussian unitary ensemble 
(GUE). Namely, it appears in the level spacing distri- 
bution at the edge of the spectrum 24|, |25( ■ The bulk 
eigenvalue density of N x N GUE matrices is described 
by the so-called Wigner semicircle with the endpoints at 
±y/2N. However, in order to capture the fine structure of 
the edge of the spectrum one has to magnify it by choos- 



ing the scaling variable 



2-1/2^-1/6^ Then the 



probability that exactly n eigenvalues lie in the interval 
(s, oo) is given by the expression [3, 



E(n,s) 



(-l) n d" 
n\ dz' 



■ det(l - zK) 
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We can easily prove now that the eigenvalue statistics 
(fT3| at the edge of the GUE spectrum is identical to that 
of the quantum front with n being the number of parti- 
cles. Indeed, the generating function of ([T3"f is obtained 
by a Fourier transform and can be written as 



exp 



-e lA — det(l-zin 
az 



(14) 
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Here, one can notice the operator of the translation which 
acts as exp(a^)/(z) = f(z + a) on any analytical func- 
tion. Translating the argument and setting z = 1 after- 
wards, one immediately recovers the previous result (1121) . 
We have thus mapped the FCS of the quantum front to 
the well studied problem of GUE edge eigenvalue statis- 
tics. One should point out, however, that while the edge 
region of the front widens with time, the edge of the ran- 
dom matrix spectrum shrinks with the matrix size and 
thus the correspondence is valid only after proper rescal- 
ing. 

As an application of the mapping, we revisit the ques- 
tion of interpretation of the edge density profile in terms 
of single particles ljj , and we also obtain the order statis- 
tics of the particles. In the random matrix language, the 
probability density of the nth largest eigenvalue for GUE 
random matrices is given by 24 1 



F(n,x) 



k=0 



dE(k,x) 
dx 



(15) 



In addition, one has the following sum rules [2£ 



X>( fc > 



J2kE(k,: 



TrK 
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with obvious probabilistic interpretations. Using the 
above sum rules, one arrives at the scaled density of 
eigenvalues as a sum of single eigenvalue densities 



J2F(n,x) = ~J2k 



k=0 



dE(k,x) 
da; 



P {x) (17) 



where 



p(x) = K(x,x) = (Ai'(.x)) 2 - xAi 2 (x) 



(18) 



is obtained from Eq. (jlip by taking the limit x — > y. 

The results ([T5|) and (fT?) have clear interpretations in 
the particle picture. Namely, the probability distribution 
of the nth eigenvalue F(n, x) is the probability distribu- 
tion of the scaled position of the nth particle. Thus the 
functions F(n, x) provide us with the order statistics of 
the rightmost particles in the front with F(l,x) being the 
Tracy- Widom distribution (24[. Furthermore, the sum of 
the probability distributions F(n, x) gives us the scaled 
density p(x) of particles in the front. Fig. [2]shows the dis- 
tributions of the first six particles obtained by a powerful 
numerical method for the evaluation of Fredholm deter- 
minants [2^1 . The sum of them, shown by the dashed 
line, deviates from p(x) only at x ~ —8. 

It is remarkable that the description of the front we 
found is consistent with a classical particle picture. One 
should also remember, however, that the position of the 
particles spreads out as i 1 / 3 in the original (unsealed) 
variables. This subdiffusive spreading has a quantum 
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FIG. 2: Rescaled density near the edge of the front for t = 100 
(diamonds) and t = 1000 (circles) . The solid (red) line shows 
the scaling function p(x) while the distribution functions of 
the nth particle for n = 1,...,6 are shown by the dotted 
(blue) lines. The dashed (green) line is the sum of the F(n, x). 



origin, namely it follows from the cubic nonlinearity of 
the dispersion around the Fermi points. 

Additional information on the front region is contained 
in the cumulants K n of the FCS. They can be obtained 
from logarithmic derivatives n n = (— id\) n lnx(A, s)| A=0 
of the generating function (fl~2"|) . Using properties of Fred- 
holm determinants, the cumulant generating function 
reads 



lnx(A,s) 



TV K k 

fc=l 



-(e l 



1) 



(19) 



where the trace is defined as 

TrK k = / dx 1 ...dx k K(x u x 2 )...K(x k ,x 1 ) (20) 

J s 

The first two cumulants, corresponding to the total num- 
ber and to the fluctuations of the particle number, have 
simple forms n\ = Tr K and n 2 — Tr K(l — K), respec- 
tively. In general, carrying out the traces in K n is diffi- 
cult and would require knowledge of the eigenvalues of 
the kernel. Having the spectrum would also give access 
to more complicated quantities of interest such as the 
entanglement entropy }10h12| 



S=-Tr[KlnK+(l-K)\n{l-K)}. (21) 

Even though there exists a differential operator commut- 
ing with K 24J, the solution of its eigenproblem is not 
known and the analytic calculation of k 2 and S remains 
a challenge. 

Numerically, we can calculate n 2 and S from the ma- 
trix representation of the FCS given in ©, using matrices 
of sufficiently large size. The results are shown in Fig. 



4 



[3] for different times, plotted against the scaling variable 
s. One can see that the convergence to the t — > oo limit 
is fast. Indeed, for t = 1000, we have a nearly perfect 
collapse onto the scaling functions K2 = Tr K (1 — K) and 
S given by Ecl (l!?Tj) and evaluated by using the methods 
described in [26|. This gives us a further check on the 
scaling function of the FCS (p"2")l . One should also note 
that both K2 and S inherits signatures of the discrete 
particle picture developed above. 




FIG. 3: Entanglement entropy S (upper symbols) and 
particle-number fluctuations K2 (lower symbols) in the edge 
region for different times as functions of the scaled coordinate 
s. The solid (red) lines are the t — > 00 scaling functions. 

We have found a direct correspondence between count- 
ing statistics of free fermions at the edge of an evolving 
quantum front and that of random matrix eigenvalues in 
GUE edge spectra. A similar correspondence between 
the equilibrium FCS of fermions in a line segment and 
the statistics of bulk eigenvalues in the GUE spectrum 
can also be recognized. Indeed, in a proper continuum 
limit, the FCS can again be written as the Fredholm de- 
terminant (|12p with the Airy kernel K(x,y) replaced by 



the sine kernel 27 



28] which is associated with the bulk 
spectra of GUE matrices. For fermions on a finite chain, 
the connection to random matrices was already pointed 
out in (29| where the Toeplitz determinant arising in the 
study of entanglement entropy was related to the eigen- 
value statistics of the circular unitary ensemble (CUE). 
Note, however, that the CUE and the GUE statistics be- 
come identical in the bulk scaling limit 25] . 

Interestingly, the equilibrium statistics also appears in 
the context of our non-equilibrium problem if one consid- 
ers a finite segment in the middle of the chain (m, n <C t) 
in the limit t — > 00. Then the matrix elements C mn (t) 
become independent of t and unitarily equivalent to the 
equilibrium correlations (l6j . as can be seen from an 
asymptotic expansion of the Bessel functions in (JT]). As 
discussed in the previous paragraph, this implies bulk 
GUE statistics for the FCS. 



It would be important to check the universality of the 
edge behaviour of quantum fronts. The simplest general- 
ization would be to start from an initial state, where the 
left (right) hand side of the chain is not completely filled 
(empty) [17j or the initial density profile is a smooth func- 
tion [21|. In such cases the staircase structure has been 
found to be essentially unchanged. Moreover, a very sim- 
ilar structure of the magnetization front was observed in 
the transverse Ising model starting from a domain wall 
initial state . One might expect that the FCS remains 
unchanged in the above cases. 

Another important question is the role of interactions. 
Does the edge structure survive if one remains in the inte- 
grable Luttinger liquid regime? How robust is it against 
integrability breaking terms? Recently, a numerical tech- 
nique was proposed [31| which might be suitable to attack 
these questions. 

Finally, we mention an intr iguing connection to the 



asymmetric exclusion process [32H34I ]. There, starting 
from a step initial condition, the distribution of parti- 
cle positions can also be calculated [35j. One observes a 
i 1 / 2 scaling for the width of the distributions at the edge 
of the front, while the t 1 ^ 3 scaling and the Tracy- Widom 
distribution emerge towards the bulk where the exclusion 
interaction between the particles becomes more impor- 
tant. Although there is no one-to-one correspondence to 
our result, nevertheless one may ask whether the quan- 
tum effects due to the curvature in the spectrum could be 
described in terms of some generalized semiclassical pic- 
ture, e.g. by introducing effective interactions between 
particles of different velocities. 
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